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The coherent transport properties of armchair graphene nanoconstrictions(GNC) are studied using 
tight-binding approach and Green's function method. We find a non-bonding state at zero Fermi 
energy which results in a zero conductance valley, when a single vacancy locates at y = 3n ± 1 of 
a perfect metallic armchair graphene nanoribbon(aGNR). However, the non-bonding state doesn't 
exist when a vacancy locates at y=3n, and the conductance behavior of lowest conducting channel 
will not be affected by the vacancy. For the square-shaped armchair GNC consisting of three metallic 
aGNR segments, resonant tunneling behavior is observed in the single channel energy region. We 
find that the presence of localized edge state locating at the zigzag boundary can affect the resonant 
tunneling severely. A simplified one dimensional model is put forward at last, which explains the 
resonant tunneling behavior of armchair GNC very well. 



I. INTRODUCTION 

Owing to many good properties, such as the ultra high 
Fermi velocity(10 6 m/s), the stability due to the sp2 hy- 
bridization and the feasibility of large-scale integration, 
graphene has been regarded as a promising candidate 
material for post-silicon electronics, and has stirred in- 
tensive studies 0, i, |, 0, i, i, 0, S II M M O ■ Many 
graphene microstructures have been studied both theo- 
retically and experimentally for future applications, such 
as the field-effect transistors, the p-n junctions, metallic- 
metallic (semiconducting) nano junctions, L-shaped, Z- 
shaped, T-shaped and cross shap ed ,junctions|l3l. 0, EH, 

m [13, tn, ii m m m, m ise 11 m 

Graphene nanoconstrictions(GNC) is one of the im- 
portant building blocks for the carbon based electronic 
circuits[H [H M, H, H H [H, H, S3] . With the com- 
bination of etching and deposition techniques, a com- 
plete turn off of electrical transport of GNC has been 
demonstrated as a function of the local gate voltage 
experimentally [31|. Furthermore, stable and rigid car- 
bon atomic chains were experimentally realized by re- 
moving carbon atoms row by row from graphene, which 
shows the sophisticated technique to control the shape of 
graphene nanodevicesjll]. Theoretical researches mainly 
focus on two types of GNCs, i.e., the type based on zigzag 
graphene nanoribbon(zGNR), and the other type based 
on armchair graphene nanoribbon(aGNR). The wedge- 
shaped GNCs based on zGNR show a gap in the trans- 
mission spectrums, and the zero conductance is related 
to the appearance of localized zero energy edge states [29|. 
Based on a zigzag edged GNC, valley filter is proposed 
which can produce a valley polarized current. Two valley 
filters in series ma y fu nction as an electrostatically con- 
trolled valley valve[30j]. As for the GNCs based on aGNR, 
performance limits of GNR field-effect transistors are in- 
vestigated by studying a GNC with two semiconduct- 
ing wide aGNRs attached to a narrow semiconducting 
one [13]. 

In this paper, we first discuss the lattice vacancy effects 
on the transport properties of metallic aGNR. We then 



study the coherent transport properties of GNCs with 
two wide metallic aGNRs attached to a narrow metal- 
lic one with two vertical zigzag boundaries (ZBs). We 
investigate the role of the ZBs of GNCs in the coher- 
ent transport by breaking the ZBs one by one. Using a 
one dimensional model for the GNCs, we are able to un- 
derstand the resonant tunneling behaviors of the GNCs. 
Throughout the paper, we use the nearest tight-binding 
approach and the recursive Green's function method to 
analyze the transmission rate and the local density of 
states(LDOS) [13, H S3, El ■ 

The paper is organized as follows. In Sec. [Til we discuss 
the lattice vacancy effects on the transport properties of 
metallic aGNRs. In Sec. IIII1 the transport properties of 
armchair GNC are studied, and the simplified model is 
put forward. In Sec. IIV1 the summary is given. 



II. SINGLE LATTICE VACANCY EFFECT ON 
METALLIC AGNR 

In the nearest tight-binding model, a perfect zGNR is 
always metallic because of the finite overlap of the two 
edge states 49]. However, aGNR can be either metallic 
(for N — 3n— 1) or semiconducting (otherwise), where N 
is the width of of aGNR counted by the number of dim- 
mer lines in the transverse direction with n an integer [lfj|. 
It has been shown that a single vacancy at the edge of 
aGNR would induce a zero conductance dip in the con- 
ductance spectrum|12]. It has also been found that an 
on-site defect at different positions of aGNR induce dif- 
ferent conductance spectrums [2(|. It is then interesting 
to investigate the transport property in the presence of 
a single vacancy at various positions of aGNR. 

The tight- binding Hamitonian of the aGNR is given by 



HaGNR = ^ Sicjci + t c\ci 



(1) 



where £j is the on-site energy of atom i, cj(cj) the cre- 
ation(annihilation) operator at site i, and t the hopping 
constant between the nearest neighbors. 
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A lattice vacancy can be described by a v ery large i 
site energy in the tight-binding model (4l|, [52] • Due 
the translational invariance of a perfect aGNR in th< 
axis, we only have to consider the effect of position 
the vacancy along the y axis. Here we choose the wic 
N = 8 as an example, see Fig. [TJ The vacancy is marl 
by a square. With the consideration of symmetry ale 
the y axis, we only have to consider the vacancy at 1 
four positions, i.e. , y = 1, 2, 3 and 4. 

The conductance spectrums are shown in Fig. [TJ 1 
solid line is the conductance spectrum for a perf 
aGNR, which shows quantized conductance platea 
The dashed line, dash dotted line, and solid line w 
solid spheres represent the cases that a single vacancy 
sides at y = 1, y = 2, and y = 4. All the three lines sh 
a conductance valley with a zero conductance at E = 
The solid line with open circles shows the conductai 
spectrum when a single vacancy locates at y = 3. We 
that the conductance spectrum maintains perfectly 1 
first quantized conduction plateau. The presence of 1 
vacancy doesn't affect the lowest conducting channel 
metallic aGNR at all. 



(a) Vacancy at y=1 




(b) Vacancy at y=2 




(c) Vacancy at y=3 




FIG. 1: (Color online). Conductance spectrums for a perfect 
aGNR(solid line) with width 8, and for the case that a single 
vacancy locates at y=l(dashed line), y=2(dash dotted line), 
y=3(solid line with open circle), and y=4(solid line with solid 
sphere). The conductance is in units of Go = 2e 2 /h, where h 
is the Plank constant. 



To understand the peculiar effects of single lattice va- 
cancy on the metallic aGNR, we examine the microscopic 
distribution of LDOS of the ribbon. Figure [Ha) shows 
the LDOS distribution at E = when a single vacancy 
locates at y = 1, as marked by the square, where the 
magnitude of the LDOS is denoted by the radius of the 
solid circle. We can see that for the sublattice to which 
the vacancy belongs (say A sublattice), all the atoms have 
zero LDOS. Finite LDOS appear only on the other sub- 
lattice^ sublattice). ft is this non-bonding state that 
gives rise to the zero conductance at E = 0, as shown in 




FIG. 2: (Color online). The LDS distribution of perfect 
aGNR (a), and the aGNR with a vacancy at y = 1,E = Ot 
(b), y = 2, E = Ot (c), y = 3, E = Ot (d), with the width N=8 



Fig. [TJ When the energy deviates from zero, LDOS on A 
sublattice becomes finite, which opens possible hopping 
paths for electrons to cross the vacancy region, and thus 
provides finite conductance. Figure^b) shows the LDOS 
distribution when the vacancy locates at y = 2. We can 
see that the non-bonding state also appears. When the 
vacancy locates at y = 4, the situation is similar with 
that of y — 1, and y = 2. However, for a vacancy at 
y = 3, we see in Fig. [2jc) that the LDOS distribution 
is uniform at y — 3n ± I, and zero at y = 3n. The 
pattern is the same throughout the energy region of the 
first conductance channel, which explains why the lowest 
conductance plateau is unaffected by the presence of a 
vacancy at y — 3 shown in Fig. [TJ 

We first notice that the response of the system to a 
vacancy at y — 3n can be explained by the wave function 
of a perfect aGNR [UBS!, 



$fc (y) 
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\/N + 1 V se 
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(2) 



with s — ±1 for electrons and holes, #(k) the phase differ- 
ence between two sublattices. The wave function exhibits 
nodes at y — 3n, which makes these sites insensitive to 
introduction of vacancy. 
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In order to understand the phenomenon that a single 
vacancy reduces the LDOS at all the atoms of the same 
sublattice to zero at E = 0, we go back to the Shrodinger 
equation for the tight-binding model. The wave functions 
at the three nearest neighbors of the atom i should satisfy 
the relation 



At zero energy, this implies 



= 0. 



(3) 



(4) 



Without lost of generality, we take aGNR with width 



X X+l 




FIG. 3: (Color online). 
aGNR. 



A single vacancy at the edge of the 



N = 8 with a vacancy on the B sublattice, as shown by 
Fig. [3J Now we set a vacancy at the edge as marked 
by the square in Fig. [3l which renders $1 = 0. From 
the Shrodinger equation we can immediately get $2 = 0. 
Then from $3 = 0, due to the node of wave function of 
aGNR, and $2 = 0, we can deduce $4 = 0. In a similar 
way we can show that all the atoms of sublattice B have 
zero LDOS. It is easy to see that the same discussion 
applies for the vacancy at any positions of y = 3n ± 1 of 
the aGNR. 

The appearance of such non-bonding state is closely 
related with the topological structure of graphene. A 
single vacancy cannot induce such non-bonding state for 
the square lattice, as shown by Fig. [H where the triangle 
denotes the vacancy. 
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FIG. 4: (Color online). LDOS distribution of square lattice 
ribbon without (a) and (b) with a single vacancy. The triangle 
denotes the vacancy. 




Lc 



FIG. 5: Schematic diagram of the square shaped GNC. N = 
11 and N c — 5 denote the widths of the wide aGNRs and 
narrow aGNR respectively. L c — 4 is the length of central 
narrow aGNR in units of a. 



will focus on the energy regime AE2 around the zero en- 
ergy, in which both the wide and narrow aGNRs exhibit 
the single conduction channel. As in Fig. [7] we observe 
oscillating conductance with the peak value reaching Go, 
and a larger L c induces more rapid oscillations. Around 
E = 0, the conduction is suppressed to zero for all L c > 1 



III. TRANSPORT PROPERTIES OF 
ARMCHAIR GNC 

We now turn our attention to the armchair GNC shown 
in Fig. [5j with widths of the wide metallic aGNRs and 
the narrow one denoted by TV and N c , and the length of 
the central segment L c in units a. Two vertical aligned 
zigzag boundaries are emphasized by the bold lines with 
the notation 61, 62. We first study the conductance as 
a function of L c . The band structures of aGNR with 
width 5 and 11 are shown in Fig. [5J In this paper, we 



To understand the peculiar behaviors of the conduc- 
tance of the GNC, we study the LDOS distribution. As 
in Fig. [SJa), there is a localized edge state at each ZB. 
We have also investigated ZB's with different heights (4, 
5, 6, 7, 8, 9 for GNC with wider aGNR N = 23) and 
found that a ZB with three successive sites in a sequence 
is enough to induce localized edge state, consistent with 
previous works [1, HJ, The localized edge state is 
suppressed by the introduction of a single vacancy at the 
ZB bl and b2, as shown in Figs.JSJb) and (c). 

Next we study how the localized edge states affect the 
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N=11, Nc=5, Lc=4 



(a) 







FIG. 7: (Color online). Conductance spectrums of square 
shaped GNC for N = 11 and iV c = 5, and L c equals to 4(a), 
8(b), and 12(c). 



transport behaviors by comparing the conductance spec- 
trums of the three GNCs shown in Fig. [8] As shown in 
Figs. US (a) and (b), for the GNCs with one or two ZBs, the 
conductance is suppressed to zero near the zero Fermi en- 
ergy. However, for a GNC without ZB, the conductance 
remains finite at E = 0, see Fig. [H] (c). Therefore, the lo- 
calized state at the ZB is responsible for the suppression 
of conductance to zero at E = 0. 

Since we are considering the energy region where both 
the wide aGNRs and the narrow one permit only one 
conducting channel, they can be modeled by a sin- 
gle nanowire. A ZB which induces localized state can 
be modeled by a quantum dot (QD) coupled to the 
nanowire. The model system is shown in Fig. 1101 Here we 
would like to discuss the simpler case first, i.e., the model 
with a single QD. For this system, Orellana et aZ.[45[ have 




FIG. 8: (Color online). LDOS distribution for a GNC with 
two ZBs (a), one ZB (b), and no ZB (c). The radius of the 
solid circle (red) stands for the magnitude of LDS. The blank 
ring(blue) in panel (b) and (c) means the vacancy. N = 11, 
Nc = 5, Lc = 4. 



already solved the transmission coefficient by setting the 
wave function as 



4, 



fe- 



re 

ikx 



x < 
x > 



(5) 



and matching the wave functions from both the left and 
right sides to that coupled by the QD. The conductance 
can be derived analytically as 



2e 2 
h 



2e 2 



(^i) 2 +(™r 



(6) 



where e = 2v cos ka is the Fermi energy of the quantum 
wire, Ei is the eigen energy of the side coupled by QD, 
v is the hopping constant both within the nanowire and 
between the nanowire and QD. From Eq. © one can 
find that an antiresonance appears at the energy s = e±. 
The localized edge state at zero Fermi energy [7J implies 
£i = 0, and explains why a single ZB of GNC induces a 
zero conductance at the zero Fermi energy. Figure. E^a) 
shows the transmission spectrum of GNC with one ZB, 
which is well described by the simple model especially in 
the small k region. 

Similarly, the GNC with two ZBs, as shown in Fig. [8] 
(a), can be modeled by two QDs coupled to a nanowire, 
as shown by Fig. [TU] (b). Now we extend the approach 
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Since c diverges at e — E\ = 0, 



where c — - ( e 

w y e— ei 

one has i = 0, and thus the transmission probability T = 
tt* = 0. Figure [9jb) shows the transmission spectrum 
of the model of two QDs. We find that the model can 
describe the GNC with two ZBs very well. 

In the absence of ZB in the GNCs, there will be no lo- 
calized state, and thus no antiresonance, which results in 
a normal oscillation of conductance without zero conduc- 
tance valley, as shown by Fig.[5](c). We also investigated 
GNCs with armchair boundaries, and found no localized 
edge states, as shown in Fig. [TTJ associated with finite 




FIG. 9: (Color online). Conductance spectrums for GN< 
with (a) one ZB, (b) two ZBs, and (c) no ZB. The dash 
curves are for the nanowire and QD model. Energy is nc 
malized by Asc, the half-width of the single channel regii 
(AE 2 for GNC and 2t for QD model). 





FIG. 11: (Color online). LDOS of GNC with armchair bound- 
aries, with N = 17, N c = 5 and L c = 4. 
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FIG. 10: (Color online). Schematic diagram of the model in 
which a quantum wire couples to a single QD (a) and two 
QDs (b). £i is the eigenenergy of the QD(s). 

by Orellana et a/.[45[ to the model with two QDs. The 
wave function should be 

{ e ikx _|_ re -ikx ^ x <0 
t'e ikx + r 'e~ lkx , < x < L (7) 
te lk< - x - LS> , L < x. 

Solving the wave function matching equations, we obtain 
the transmission coefficients 

2^ sin 2 k 

t = sinM _ (g) 

/ sinfc \ 2 _ ( e jk _ c , sinfc(L-l) \ 
V sin kL J \ sin kL J 



IV. SUMMARY 

Using tight-binding approach and Green's function 
method, we show that for a metallic aGNR with a sin- 
gle vacancy at y = 3n ± 1, the zero conductance valley 
arises from a non-bonding state at the zero Fermi energy. 
Next, for the square-shaped armchair GNC consisting of 
three metallic aGNR segments, we find resonant tunnel- 
ing behaviors in the energy regime with single conduc- 
tion channel. It is shown that the localized edge state 
at the ZB affects the transport properties severely. A 
one dimensional model with one or two QDs coupled to 
a nanowire is put forward which can explain the resonant 
tunneling behavior of aGNC very well. 
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